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SUMMARY 

Structural technology of laminated filamentary-composite stiffened-panel structures 
under combined inplane and lateral loadings is discussed. Attention is focused on (1) 
methods for analyzing the behavior of these structures under load and for determining 
appropriate structur.,1 proportions for weight-efficient configurations, and (2) effects of 
impact damage and geometric imperfections on structural performance. Recent improvements 
in buckling analysis involving combined inplane compression and shear loadings and trans- 
verse shear deformations are presented. A computer code is described for proportioning or 
sizing laminate layers and cross-sectional dimensions, and the code is used to develop 
structural efficiency data for a variety of configurations, loading conditions, and 
constraint conditions. Experimental data on Luckling of panels under inplane compression 
is presented to validate the analysis and sizing methods and to illustrate structural per- 
formance and efficiency obtained from representative structures. Experimental results 
show that strength of panels under inplane compression can be degraded by low-velocity 
impact damage. Mechanisms of impact-damage initiation and propagation are described. 
Finally, data are presented that indicates the matrix is a significant factor influencing 
tolerance to impact damage. 


INTRODUCTION 

To take advantage of the mass-saving potential of advanced composite materials, pro- 
cedures must be developed that provide reliable and efficient structural designs. In 
addition to satisfying traditional design requirements typical of metal structures, 
composite structural designs must also account for characteristics and failure modes uni- 
que to composite materials. Once these features are understood, analytical methods can be 
developed to predict the behavior of structu; ‘ components made of composites, and appro- 
priate criteria can be imposed to insure th. efficient composite structures meet all 
design requirements. 

One research focus at the NASA Langley Research Center has been the development of 
analysis and sizing methods for composite structural panels required to carry compression 
and combined loads. A state-of-the-art review of these studies was presented in 1975 
which included an investigation of panel buckling as well as the effects of low-velocity 
impact damage on composite sandwich panels (ref. i). During the past three years, consi- 
derable progress has been made in understanding ti.e structural behavior of stiffened 
composite panels and several failure mechanisms that affect the performance of compression 
panels subjected to low-velocity impact damage have been identified. The graphite-epoxy 


) 


S' a? 


materials from which the panels were fabricated are commercially available 450K cure 
temperature systems which were processed in an autoclave following the manufacturer's 
recommended procedures. The present paper summarises Langley research activities conduct- 
ed since 1975 on stiffened composite panels. Advances in analysis and sizing procedures 
for stiffened compression panels are discussed, and experiments conducted to verify these 
analysis and siting procedures are described. The effect of low-velocity impact damage on 
the strength of compression panels is also presented. 
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SYMBOLS 

Planform area of stiffened panel. 

Width. 

Young's modulus of composite material in fiber direction. 
Longitudinal extensional stiffness of panel. 

Amplitude of overall bow at panel midlength. 

Transverse shear modulus. 

Shear stiffness of panel. 

Panel length. 

Stress resultants. 

Load index. 

Lateral pressure loading on panel. 

Thickness . 

Mass of stiffened panel. 

Mass index. 

Amplitude of eccentricity at panel midwidth. 

Strain at buckling. 


STIFFENED PANEL ANALYSIS AND SIZING 

The complexities of laminated composite panels have led to development of sophisti- 
cated analysis and design procedures compared to those formerly used for metal panels. 
Panels must be designed to carry combined loads and structural efficiency and manufactur- 
ing considerations require that a variety of structural configurations be considered 
(see figure 1). The VIPASA computer code (ref. 2) is an example of a sophisticated 
analysis program for the rapid and accurate buckling solution of stiffened composite 
panels subjected to compression loading. VIPASA, however, gives conservative solutions 
for the general buckling mode of a rectangular panel subjected to shear loading and does 
not include the effects of through-the-thickness transverse shear deformations. Approaches 
for improving the solution for both of these problems, and typical examples are described 
herein. In addition, a sizing code called PASCO (Panel Analysis and Sizing Code) which 
uses VIPASA to perform buckling analyses is described, and the PASCO code is used to pro- 
vide structural efficiency results for important generic classes of problems. 
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Figure 1.- Stiffened panel subjected to combined loads 
and representative stiffener configurations. 



Figure 2.- Comparison of predicted buckling loads from various analyses for 
blade-stiffened panel subjected to combined longitudinal compression and shear loadings. 
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Stiffened Panel Stability Analysis 


Combined shear and compression loads .- The VIPASA buckling analysis (ref. 2) provides 
an exact solution to the classical thin-plate equations which are based on the Kirchhoff- 
Love hypothesis. The buckling solution is obtained by assuming a sinusoidal buckle pattern 
in the stiffener direction. A sinusoidal buckle pattern assumption for orthotropic panels 
loaded by inplane longitudinal and transverse loads permits simple support boundary condi- 
tions to be satisfied on edges normal to th< stiffeners while the boundary condition on 
edges parallel to the stiffeners may be arbitrarily specified. The presence of shear or 
anisotropy causes skewed node lines which do not conform to the panel rectangular bound- 
aries. In these situations, the resulting theoretical solutions are lower than that given 
by the solution from a two-dimensional analysts in which the desired boundary conditions 
are accurately modeled. The error introduced by anisotropy for most practical configura- 
tions is small, but the buckling strain for the general mode (m ■ 1) of a panel loaded in 
shear can be substantially in error. Buckling solutions for local modes in which more than 
one buckle forms along the length, however, are accurately predicted by VIPASA. 

Because of the complexity and computer expense of a full two-dimensional analysis, an 
approximate method based on a combination of VIPASA and orthotropic plate theory has been 
developed for panels subjected to combined shear and compression loading. Typical inter- 
action curves for shear and compression are shown in figure 2 for a 76.2-cm-square, blade- 
stiffened composite panel having six stiffeners. Materials properties used in the analysis 
are presenr 'd in Table I (material A-tape) and cross sectional dimensions are defined in 
Table II (design B). The desired boundary conditions are simple support on all four 
edges. In the VIPASA analysis, however, boundary conditions can be specified only on the 
two edges parallel to the stiffeners. The results of this analysis are presented as the 
solid curve in figure 2. VIPASA was also used to obtain the buckling load of the panel 
simply supported along the side edges assuming smeared otthotropic stiffness properties 
(lower dash curve). Differences in these two results are attributed to the inadequacy of 
representing the stiffened panel by average orthotropic stiffnesses. The same orthotropic 
panel can be analyzed by VIPASA with the edges perpendicular to the stiffeners simply 
supported. This result (the upper dashed curve) shows a large increase in buckling load 
reflecting the fact that it is more important to model correctly the boundary conditions 
on the edge normal to the stiffening. It is postulated that any difference in the lower 
two curves would be similar to the difference between the upper curve for the VIPASA 
analyses using orthotropic stiffnesses and an exact analysis of the detailed cross section 
with ends simply supported. Thus, the ratio of the two lower results is applied to the 
upper curve to give the curve labeled Adjusted VIPASA. 

The accuracy of the approach for predicting the buckling of rectangular panels loaded 
in shear is indicated by a few results obtained from the STAGS computer program (ref. 3) 
which treated the detailed panel simply supported on all four edges. Either of the two 
upper curves provide a conservative but reasonably accurate estimate of the correct result 
obtained from the two dimensional STAGS analysis. Not satisfying boundary conditions on 
the side edges of stiffened panels leads to little error because the node lines tend to 
align approximately with the stiffeners and several buckles form across the width. For 
certain configurations where local stiffener deformation is pronounced, the adjusted VIPASA 
result can be significantly lower than the orthotropic plate results with simply supported 
edges. The conservative approach of choosing the lower of these two approximate solutions 
is used herein. Additional studies have been made with the approximate analysis 
that show it to be in reasonable agreement with more accurate analyses whenever the buckle 
length transverse to the stiffeners is greater than twice the stiffener spacing. In the 
example shown in figure 2, the buckle length transverse to the stiffeners was approximately 
3 stiffener spacings. Design studies on panels loaded by shear and biaxial compression 
have shown that using the conservative VIPASA solution leads to only a few percent mass 
penalty for most practical load combinations. If shear is the dominant loading, the 
approximate analysis can be used to take advantage of the greater predicted load carrying 
capability if proper caution is used in evaluating the results based on the calculated 
buckling mode shapes. The advantages of Che approximate analysis relative to a two- 
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TABLE I.- GRAPHITE-EPOXY LAMINA PROPERTIES 


MATERIAL A MATERIAL B 



Tape 

Fabric 

Tape 

Modulus in fiber direction, GPa 

131.0 

62.7 

110.0 

Modulus normal to fibers, GPa 

13.0 

62.7 

- 

Shear modulus, GPa 

6.4 

6.4 

- 

Major Poisson's ratio 

0.38 

0.10 

- 

3 

Density, Mg/m 

1.52 

1.52 

1.60 

Thickness, mm/ply 

0.14 

0.36 

.17 


dimensional analysis are (1) less effort is required in data preparation, (2) a solution 
convergence check is not required and (3) three orders of magnitude less computer time is 
required. 

Transverse shear eff ects . - Composite panels with open-section stiffeners such as 
blade stiffeners have been shown (ref. 4) to require a higher mass to carry a specified 
compression load than stiffened panels with closed-section stiffeners such as hat 
stiffeners. One open-section stiffener configuration which has been shown to have imp- "^ved 
structural efficiency is the sandwich-blade (ref. 5) in which the web is of sandwich on- 
struction and connects the panel skin to a cap composed mostly of 0° oriented plies. 

However, some of the advantage of this configuration is offset by the low transverse shear 
stiffness of the honeycomb core. 

The VIPASA analysis does not include transverse shear deformations, although finite 
element codes such as NASTRAN (ref. 6) have plate elements which do include transverse 
shear deformations in the buckling formulation. However, an analysis more rapid than 
NASTRAN is desired that can be incorporated in a structural sizing code. One such analysis 
has been developed and is presented in reference 5. 

An illustration taken from reference 5 showing the influence of transverse shear on 

buckling is presented in figure 3. The buckling strain is presented as a function of the 

buckle length for designs having a honeycomb core with transverse shear modulus values of 
0.26 GPa, 1.21 GPa, and for an infinite transverse shear modulus. The solution with infinite 
transverse shear modulus was obtained using the BUCLASP2 computer code (ref. 7). The 
dimensions of the cross section and the number of plies in each element are presented in 

Table II as design A. For the results shown in figure 3, three buckling modes occur 

including local buckling of the skin, twisting of the stiffener, and wide-column buckling 
of the panel. In the range of buckle length for which twisting is the buckling mode 
(10 to 70 cm), the buckling strain is strongly dependent on the transverse shear stiffness 
of the core. For a 30-cm-long buckle length, for example, the 0.26 GPa shear modulus 
reduces the buckling strain by 40 percent relative to the solution for infinite shear 
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Figure 3.- Effect of transverse shear on the buckling strain 
of a sandwich-blade stiffened panel. 
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Figure 4.- PASCO panel optimization capabilities. 
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stiffness, and the 1.21 GPa shear modulus reduces the buckling strain by IS percent. The 
NASTRAN solutions and experimental results correlate closely with the reference 5 approxi- 
mate analysis solutions. 


Stiffened Panel Sizing Procedure - PASCO 

The VIPASA buckling analysis described in the preceding section has been incorporated 
in a computerized sizing program denoted PASCO (Panel Analysis and Sizing Code). Some 
important capabilities of PASCO are indicated in figure 4. The code can be used to size 
stiffened panels having an arbitrary configuration subjected to any combination of inplane 
loadings (tension, compression and shear) and lateral pressure. The panel cross section 
is modeled as an assembly of flat plate elements in which each of the plate elements is a 
Oalaaced symmetric laminate with an arbitrary number of layers. The panel cross section 
and loading are assumed to be uniform in the direction of the stiffeners. Stresses 
associated with a bow-type initial imperfection or lateral pressure are accounted for 
using a beam-column approach. 

W/A 

During sizing, the merit function is the mass index — , the mass per unit area o'" 

Ij 

the panel divided by the panel length. The sizing variables are the individual plate 
widths, ply thicknesses, and ply orientations. Inequality constraints can be placed on 
buckling loads, lamina stresses or strains, and overall panel stiffnesses. The objective 
of the procedure is to determine the values of the sizing variables which minimize the 
mass index and satisfy the prescribed set of inequality constraints. Additional discussion 
of PASCO is presented in reference 8. 



Figure 5.- Structural efficiency of various stiffened panel configurations. 
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PASCO APPLICATION EXAMPLES 

In the two studies presented below denoted "longitudinal compression" and "combined 
loads", the panels were assumed to be geometrically perfect and the only imposed constraint 
was buckling. The increases in mass necessary to meet the additional requirements imposed 
by a geometric imperfection and by prescribed extensional and shear stiffnesses are presented 
in the section entitled "effect of overall bow and stiffness requirements." 

Longitud inal co mpres sion.- Structural efficiency data for graphite-epoxy panels with 
several cross sectional configurations sized to carry only longitudinal compressive loads, 
Njj, are presented in figure 5. Also shown for comparison is a curve for the minimum mass 
of hat-stiffened panels constructed of aluminum which was originally presented in ref. 9. 

The least efficient graphite-epoxy panels are the blade- and 1-stiffened panels. These 
two panel configurations have roughly the same efficiency, and are about 40 percent lighter 
than aluminum hat-stiffened panels. The curve for panels having a honeycomb sandwich-blade 
are approximately 20 percent lighter than panels with a solid blade. Because the trans- 
verse shear flexibility required for these panels cannot be modeled with the present version 
of PASCO, sandwich-blade stiffened panels were sized using the analysis of reference 5 
which includes the transverse shear effects. Finally, the lightest panels shown are the 
graphite-epoxy hat-stiffened panels (ref. 8). These panels weigh 60 percent less than 
aluminum hat-stiffened panels. 

C ombined loads .- Blade-stiffened panels having the configuration and variables 
described in reference 8 were sized for pure longitudinal compression and for selected 
ratios of biaxial compression and shear. The panels were square (76.2-cm-length sides) 
and were sized to have an integer number of stiffeners. The results are presented in 
figure 6. The optimum number of stiffeners varies with loading from 16 stiffeners for 
lightly-loaded panels (0.07 to 0.10 MPa) to eight stiffeners for heavily-loaded panels 
(5 to 7 MPa). 

E ffect of overall bow and stiffness requirements .- One type of geometric initial 
imperfection that is common in panels and that can be accounted for with PASCO is the over- 
all bow-type geometric imperfection. The first-order effect of the bow imperfection is 
assumed to be the additional stresses produced by bending. In I'ASCO, these bending 
stresses are calculated using a beam-column approach and are added to the stresses produced 
by the inplane loading. The resultant stresses are used to calculate the buckling ’oads 
and are also examined for maximum stress and strain limitations. 

Panels also ..onanonly have requirements which can be met by imposing limits on stiff- 
ness. Several smeared orthotropic stiffnesses are calculated in PASCO and can be used as 
inequality constraints during sizing. Studies were carrii r out to determine the effects 
of a bow imperfection and of extensional and shear stiffness requirements on the mass of 
minimum-mass blade-stiffened panels. The results are presented in figure 7. The panels 
have the same conf igu''ation as those presented in figure 6 except that all panels have 
eight stiffeners and no limitation is placed on panel width. The lowest curve is for 
panels having neither a bow nor a stiffness requirement. The next higher curve is for 
panels designed for an overall bow of e/L = +0.003, where e is the eccentricity at panel 
midlength, L is the panel length, and the + sign means that panels carry the load 
whether the bow is positive or negative. The highest curve shows the effect of adding 
shear and extensional stiffness requirements typical of transport aircraft wing panels. 

The stiffness requirements are taken from reference 1 and are given in Table III. These 
panels are also assumed to have an overall bow of e/L = +0.003. The results indicate that 
designing for an overall bow of e/L = +0.003 causes a mass penalty of about 15 percent 
and that adding a representative stiffness requirement results in an additional penalty 
of 18 to 50 percent. 
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Y- MPa 

Figure 6.- Structural efficient^ of graphite-epoxy blade-stiffened panels 
designed for specified ratios of combined inplane loads. 

Panels are square, 76 cm on a side. 



Figure 7.- Effect of selected constraints on the structural efficiency 
of graphite-epoxy blade-stiffened panels. 
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TABLE III.- MINIMUM STIFFNESS REQUIREMENTS REPRESENTATIVE 
OF COMMERCIAL AIRCRAFT WING PANELS. 


Loading Intensity Longitudinal Stiffness Shear Stiffness 


N /L® 

X ■ 

ET, 

GT, 

MPa 

MN/m 

MN/m 

0.689 

386 

79 

1.720 

508 

105 

3.440 

613 

133 

5.510 

701 

165 


^L = 76.2 era 


EXPERIMENTAL BUCKLING AND STRUCTURAL EFFICIENCY STUDIES 


Theory-Experiment Buckling Correlation 

It is important to understand the accuracy with which the buckling load of a compo- 
site panel can be predicted by theory. Discrepancies as high as 40 percent were reported 
in reference 9. The explanations given for these discrepancies were variations in thick- 
ness, variations in material properties from the nominal values used in analysis, and 
anisotropic and residual thermal stress effects not accounted for in the analysis. 

Considerable effort is required to conduct an analysis which is sufficiently accurate to 
permit valid comparisons between experiment and theory. To illustrate some of the impor- 
tant factors which must be considered, a detailed analysis was performed for a blade- 
stiffened panel loaded to buckling. The panel constructed of graphite-epoxy tape 
(material A) using design F (Table II) was •hree stiffeners wide and 76.2 cm long. 

The specimen buckled at a load of 1.76 MN/ and axial strain of 0.0037. The corres- 
ponding analytical prediction for the panel based on nominal properties and dimensions was 
a buckling load of 1.52 MN/m and axial strain of 0.0032. The effects that deviations from 
nominal values have in causing differences between experiment and theory are shown in 
figure 8. All barred quantities in the figure are associated with nominal material proper- 
ties and dimensions assumed in the design analysis. The moire-fringe pattern photograph 
(fig. 8(a)) shows the specimen buckled into four half waves along the length which is in 
agreement with the theorttical prediction (the fringe pattern for one-half wave is missing 
in the photograph due to the limited size of the grid). 

Two fundamental quantities can be measured experimentally - the imposed strain and the 
total load. To evaluate a buckling theory, comparing experimental and theoretical buck- 
ling strains is useful because the buckling strain is essentially constant with respect 
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to inplane elastic properties (see fig. 8(b)). In this figure only the dominant material 
modulus was varied while all other moduli properties were held constant at the values 
shown in Table I. The buckling strain, however, is sensitive to ply thickness variations 
as shown in the buckling strain versus skin thickness graph of figure 8(b). A thorough 
survey established that the specimen blade ply>thicknesses were essentially nominal but that 
the skin ply-thickness was 95 percent of nominal. The skin thickness variation accounts 
for a four percent reduction in the theoretical buckling strain. Transverse shear defor- 
mations of the thick blade-stiffener are responsible for further reductions in the buck- 
ling strain as shown in figure 8(c). The value of Che transverse shear modulus is not 
readily available, but considering the properties of the matrix and reported values for 
inplane shear modulus, it is probably within the range of from 30 to 60 shown on 

Che plot. Transverse shear, therefore, may account for a five Co eight percent reduction 
in Che buckling strain. The cumulative effect on buckling strain of material property 
and thickness variations from nominal compared to the experistenCal results is shown in 
figure 8(d). Correlation of theory and experiment is good when all of these factors are 
considered collectively, whereas looking at a comparison based on nominal properties would 
indicate a 13 percent discrepancy. 

A fundamental quantity of interest to the structural analyst is the structure buckling 
load. If Che material properties are not accurately known, the actual buckling Ijad may 
not agree with theory even Chough Che strain is in agreement. The buckling load is a 
direct function of Che material stiffness properties as illustrated in figure 8(e). The 
upper curve is for the nominal skin thickness and classical theory, and the lower curve is 
for a reduced skin ply thickness of t/f ° 0.93 and a correction for transverse shear 
deformation corresponding to “ 50. 

An indication of the accuracy of the assumed material properties (Table I) can be 
obtained from comparison of the theoretical and experimental axial stiffnesses shown in 
figure 8(f). The two theoretical curves are for the nominal skin thickness and a reduced 
skin ply thickness of t/t = 0.95 corresponding to Che test panel. The experimental stiff- 
ness (ET/ET = 0.99) was obtained from the applied load and the average of_several longitu- 
dinal strain gage readings. Use of the nominal material properties (En/En = 1) gives 
excellent correlation between experiment and the lower curve for the test panel 
(t/t = 0.95). Small coupon tests of samples cut from the panel also confirmed the selec- 
tion of material properties. Use of nominal material properties also permit close corre- 
lation of theoretical and experimental results for the buckling load (see the lower curve 
on figure 8(e)). Since the nominal material properties are reasonably accurate for the 
test panel, the discrepancy between the experiment and theory for buckling load based on 
nominal properties and classical theory is the same as the discrepancy for strain; namely, 
13 percent. These remits show that excellent agreement between theory and experiment 
can be obtained for both strain and load if accurate material and geometric properties and 
sufficiently precise theory are used. 

Other factors which can affect comparison of theory and experiment are boundary 
conditions and geometric imperfections. For the local buckling mode shown in figure 8, 
loaded-end boundary conditions effects were small as indicated by a finite element analysis 
on a similar configuration presented in reference 5. The effect of an overall bow imper- 
fection were also small since the bow imperfection was measured to be small (e/L = 0.001). 
In addition, a bow-type imperfection has a greater effect on the Euler buckling load which 
was significantly higher than the twisting buckling mode exhibited by this test panel. 


Panels Designed to Meet Typical Commercial- 
Aircraft Wing-Stiffness Requirements 

Experiments were conducted on graphite-epoxy panels which were designed to meet not 
only bucxling and strength constraints but also to have extensional and shear stiffness 
properties typical of coomercial aircraft wing structures (See Table III). Five panels 
(see figure 9) were tested including three blade- and tv-o hat-stiffened panels. 
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Figure 9.- Hat- and blade-stiffened panels designed 
commercial airerr^t wing stiffness requirements 
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Figure 10.- Structural efficiency of compression panels designed 
to meet commercial aircraft wing stiffness requirements. 
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Cross sectional -Jimcnsions along with the number anJ orientation of plies for characteris- 
tic elements of the cross section are presented in Table II (designs C, F, G, H and 1). 
Experimental data for the panels are presented in Table IV. The specimens were 1.47 m 
long which when flat-end tested yielded an effective simple support length of 76.2 cm. 


The experimental results for these specimens are plotted on the structural efficiency 
graph of figure 10. For comparison, the structural efficiency curve for a graphite-epoxy 
blade-stiffened panel conf igur.Ttion designed to meet bow imperfection and wing-stiffness 
requirements (reproduced from figure 7) as well as structural efficiency data for typical 
commercial aircraft aluminum wing panels (taken from reference 1) are presented. Stiffness 
is such a strong design factor for the load range shown that the theoretical curve for a 
graphite-epoxy hat-stiffened panel configuration would differ only slightly from that shown 
for a blade-stiffened configuration. The mass savings for the graphite-epoxy panels com- 
pared to the aluminum wing designs range from over 50 percent for the lightly-loaded blade 
designs to approximately 30 percent for the heavily-loaded designs. The test data for the 
light ly- loaded panels fall on or below the theoretical curve while the test data for the 
heavily-loaded panels plot above the curve. Lightly- loaded panels Cl and FI were loaded 
until they buckled, which occurred without failure. The heavily-loaded oanels Gl, HI, and 
11 were not buckling critical at the design load, and the tests were terminated at or near 
the design load to permit the specimens to be used in damage tolerance studies. The arrows 
drawn on figure 10 indicate the panels have additional load capability. 





TABLE IV.- EXPERIMENTAL DATA FOR 1.47-m-LONG PANELS DESIGNED TO MEET COMMERCIAL AIRCRAFT 
REPRESENTATIVE EXTENSIONAL AND SHEAR STIFFNESS REQUIREMENTS. TEST CONDITIONS 
SIMULATE 76.2 cm SIMPLE-SUPPORT LENGTH. 





IMPACT-DAMAGE TOLERANCE OF COIPOSITE STRUCTURES 
LOADED IN COMPRESSION 


Structural efficiency studies, such as those presented in the previous section, have 
demonstrated the mass-saving potential of high-modulus filamentary material for compression 
panels designed to meet buckling and stiffness requirements. The usefulness of composite 
materials, however, also depends on their ability to carry loads when subjected to damage 
sustained during the life of the cosiponent. The sensitivity of compression loaded struc- 
tur..'s to low-velocity impact damage was studied using flat-plate laminates as well as 
sti ' fened panels. 


Flat Laminates Damaged by Low-Velocity Impact 

Test configuration and apparatus .- The two test configurations shown in figure 11 
were used to study the effect of low-velocity impact on flat laminates. The single-bay 
specimen shown on the left is approximately 24.8 cm long and 11.4 cm wide while the three- 
bay specimen shown on the right is approximately 24.8 cm long and 38.1 cm wide. The loaded 
ends of the specimen were ground flat and loaded in fixtures that approximated clamped-end 
cor.'litions. Edge and interior restraints approximated simple-support boundary conditions. 
Thi specimen size and edge-support conditions were selected to permit the plate to have a 
I strain at buckling and therefore permit damage-initiated failure to occur at strains 
lov ;r than the buckling strain. Tiie projectile used was a 1.27-cm-diameter aluminum 
sph<:re that was propelled to impact normal to the specimen at velocities ranging from 30 
to 150 m/s. A moire-fringe technique was employed to monitor the panel lateral displace- 
ment response during loading. The front surface of the specimens was painted white to 
enhance the contract of the moire fringes. Specimens were instrumented with at least four 
stri'in gages. 

Physical characteristics of impact damage .- Specimens were impacted in the single-bay 
test frame and removed to observe the physical characteristics of damage. Following 
observation of the front and back surface damage, the specimens were ultrasonically 
inspected. Several were cross-sectioned through the impacted region and inspected micro- 
scopically for interior damage. 

Ultrasonic C-scan data and photomicrographs of the cross-section for a 48-ply laminate 
ith a {±45/02 / ± 45/02/±45/0/90}2s stacking sequence are presented in figure 12. Data 
a -e fhown for specimens impacted at 58 m/s and 95 m/s. Visual inspection showed no front 
surface damage for .-ither specimen. The specimen impacted at 58 m/s did not develop back 
surface damage -le the specimen impacted at 95 m/s had an outer ply crack on the back 
surfa:e which extended approximately 3.4 cm. Ultrasonic inspection showed both specimens 
to have sustained interior damage with the greater affected area associated with the higher 
velocity (figure 12(a)). 

Low power microscopic inspection of the cross sections in the damaged region (figure 
12(b)). shows limited damage in the specimen impacted at 58 m/s and extensive damage in 
the specimen impacted at 95 m/s. High power photomicrographs of the cross sections in the 
damaged region taken wich a scanning electron microscope are shown in figure 12(c). Both 
specimens exhibit telamination cracks between plies of dissimilar orientation (i.e., 0/90, 
0/45, and +45/-45) .s well as intraply cracking ( through-the-thickness) of 0° and ±45° 
oriented plies. Failure appears to propagate by matrix fracture and may be associated with 
cohr.iive fal ' res in the matrix or adhesive failures in the fiber-matrix bond. These 
observations show that low-velocity impact can result in significant laminate interior 
damage i.. the vicinity of the impact. Several mechanisms may participate in creating the 
local < nage illustrated in figure 12. The plate deformation response following impact is 
illu 'rated in figure 13. On contact, internal stress waves are initiated within the panel 
wh’ch may be responsible for developing local damage. The V-shaped patterns observed in 
I damaged panel (see figure 13(b)) are similar to the fracture patterns created by stress 


Mqe is 

quality 


ORIGINAL 

Of POOR 





(b) LOW POWER MAGNIFICATION, (c) HIGH POWER MAGNIFICATION 


Figure 12.- Interior damage to graphite-epoxy laminates following impact 


Figure 13.- Deformation response of laminate following low-velocity impact 
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waves in homogeneous brittle plates that have been impulsively loaded on one surface. 

Based on preliminary calculations using bulk material properties, the time required for 
translations of the initial dynamic wave through the thickness would be on the order of a 
few microseconds. Once damage has been initiated within the panel, the area in the damaged 
region has a reduced local stiffness and may deform locally out-of-plane through a combina- 
tion of stress and local bending of sublaminates. The 48*ply panels impacted at 91 m/s by 
a 1 . 27- cm-diameter aluminum projectile had local deformations in the impact damage region 
that were approximately 3 cm in diameter and a lateral deformation of approximately 1 mm 
(figure 13(c)). This local out-of-plane deformation may cause damage to propagate or 
cause additional delamination and intraply cracking to develop. It is anticipated that 
this local deformation occurs much later in time, probably several orders of magnitude 
after the initial dynamic wave effects. The out-of-plane deformation due to the plate 
overall structural response was not measured in this investigation; however, preliminary 
estimates indicate that it would have occurred later in time (figure 13(d)). The overall 
structural response probably causes little additional local damage. 

Effect of impact damage on strength .- The effect of impact damage on the strength of 
a 48-ply orthotropic laminate with a {±45/02/i45/02/i45/0/90}2s stacking sequence is 
reported in reference 10. This information along with additional unpublished test data is 
presented in figure 14. In these tests, the specimens were impacted while loaded to a 
prescribed axial strain. Specimens that failed catastrophically on impact are represented 
by filled circles, and specimens that continued to carry load after impact are represented 
by open circles. The strain at failure for control specimens are shown for comparison on 
the ordinate as open circles. Control specimens either failed in the loaded end region or 
buckled and the test was terminated. The projectile kinetic energy is plotted on the 
abscissa and discrete projectile velocities are also indicated. The curve labeled "failure 
threshold strain" in figure 14 separates results of specimens that failed catastrophically 
on impact from those that did not. The trend of the data indicates the compressive strength 
of these specimens is affected severely at the higher impact energies. For an impact 
velocity of 100 m/s, the failure threshold strain was 0.0028. 

The test results for the 58 and 95 m/s impact conditions (figures 12 and 14) suggest 
the extent of interior damage is a key factor affecting the failure threshold strain. 
Reference 10 indicates a coupling exists between the applied inplane load and the local 
deformation response associated with impact. This coupling results in a larger local 
damage region for specimens damaged by impact under load than for specimens which have no 
applied load at impact. 

Additional results are reported in references 10 and 11 of the residual strength of 
specimens which did not fail on impact. Specimens impacted at zero axial load in the 50 
to 60 m/s range were found to have ultimate residual strains significantly higher than the 
failure threshold strain. Specimen impacted at zero load in the 80 to 100 m/s range, how- 
ever, were found to have ultimate residual strains only slightly higher than the failure 
threshold value. 

Data presented in figure 14 represents specimens that were both 11.4 cm (single-bay) 
and 38.1 cm (three-bay) wide. The failure strain and failure mode was unaffected by the 
specimen width. The damage propagated laterally from the point of impact to the plate 
edge for both types of specimens. Photographs of failed panels of both specimen sizes are 
presented in figure 15. 

Damage propagation mechanisms .- Three damage propagation mechanisms for compressively 
loaded laminated composites have been observed. These include: delamination, axial load- 

lateral deformation coupling, and local shear failure. These propagation mechanisms are 
illustrated in figure 16. The delamination mechanism (figure 16(a)) can be described as 
the progressive growth of local interlaminar disbonds during the application of load. The 
propagation is dependent upon the development of local out-of-plane buckling deformations 
in the delaminated region. The closer the delaminations are to the surface, the lower is 
the load required to initiate local buckling. The series of moire-fringe photographs shown 
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of impact energy on graphite-epoxy laminate failure 
Projectiles used were 1 . 27-cm-diameter aluminum spheres 


strain 


Figure 15.- Photographs of typical failures of impact-damaged 
graphite-epoxy compression panels. 
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in figure 16(a) are typical of the progressive delamination growth which occurs with applied 
load. This failure mechanism is representative of specimens which fail during residual 
strength tests following impact and has also been observed as the failure propagation 
mechanism for impact-damaged panels subjected to cyclic compression loads (ref. 10). 

The axial load-lateral deformation coupling mechanism is illustrated in figure 16(b). 
This propagation mechanism occurred in panels impacted under load at strains near the 
failure threshold value. The local lateral deformation caused by impact, interacts with 
the applied axial load to c .use lateral propagation of the locally deformed damaged region. 
The photograph presented in figure 16(b) shows the side opposite the contact region of the 
failed three-bay panel shown previously in figure 15. Cross section A-A taken through the 
center of the failed region is shown on the right along with the 1 . 27-cm-diameter projectile 
drawn to scale. The local deformation in the impact zone and the extensive dclamination in 
the laminate may be seen in the cross section. The interior supports of the test fixture 
for the panel shown in figure 16(b) restrained the panel lateral deformations and, therefore, 
arrested the propagation of the damage zone as an axial load-lateral deformation coupling 
mechanism. The failure, however, propagated through the interior support region and to the 
panel edges by a shear failure mechanism. 

The local shear failure mechanism involves short-wavelength transverse-shear failure 
such as that illustrated in figure 16(c). The photograph on the left is a cross section 
from the outer bay of the panel shown in figure 16(b). The cross section shown on the 
right of figure 16(c) is from another specimen and illustrates more clearly the local 
shear failure mode. The damage propagation mechanism observed for a typical graphite-epoxy 
specimen is a combination of delamination and local shear. If the specimen is loaded near 
the failure threshold during impact, the failure also involves the axial load-lateral 
deformation coupling mechanism. 


Effects of Holes 


Tests were conducted on specimens containing an open circular hole to permit comparison 
of the effect of impact with that resulting from a controlled damage. Detail results of 
the study are reported in reference 10. The specimens had a {±45/02/-45/02/-45/0/90}2g 
stacking sequence with hole diameters ranging from 0.16 to 3.81 cm and corresponding 
diameter-to-width ratios ranging from 0.014 co 0.33. The effect of hole diameter on fail- 
ure strain is presented in figure 17. 

The specimen with a 0. 16-cm-diameter hole failed away from the hole at an axial com- 
pression strain on the same order as control specimens. Specimens with larger hole diameters 
failed at strains significantly lower than control specimens. For example, the specimen 
with a hole diameter of 3.81 cm failed at an axial strain of 0.0038. The failure threshold 
strain for the same laminate damaged by impact at 100 m/s by a 1 . 27-cm-diameter aluminum 
projectile was 0.0028 (see figure 14). 


A photograph of a failed specimen with a 1 . 91-cm-diameter hole is shown in figure 18(a). 
At an applied strain value of approximately 95 percent of ultimate, local delamination was 
observed at the hole edge. This delamination continued to propagate laterally to the plate 
edges with increasing load. A cross section through the hole illustrates that massive 
delamination occurred in these graphite-epoxy specimens. For comparison, a photograph is 
presented in figure 18(b) of a 0.64-cm-thick 7075 aluminum plate also containing a 1. 91-cm- 
diameter hole. The 11.4 cm by 24.8 cm aluminum specimen buckled at an axial stress of 
540 MPa. This value approaches the ultimate stress for this aluminum material and indicates 
the hole had negligible effect on reducing the specimen strength. A cross section through 
the center of the aluminum specimen is shown in figure 18(b). An increase of approximately 
11 percent in the thickness was observed locally in the region adjacent to the hole. The 
ductility of aluminum permits large local deformations to occur in this high stress concen- 
tration region without causing catastrophic failure. In comparison, the graphite-epoxy 
specimen lacks ductility and responds to the local stress con^ntration by local delamina- 
tion and specimen f.iilure. 
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(a) GRAPHITE-EPOXY (b) ALUMINUM. 

Figure 18.- Photographs of graphite-epoxy and aluminum specimens 
containing a 1 .9-cra*diameter hole loaded in compression to failure. 
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Figure 17.- Effect of holes and 100 m/s impact 
on the compression strength of a graphite-epoxy l-aminate 
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Stiffened Panfils Damaged by Impact 

The strength of structurally-ef f it ieat hat-stiffened panels damaged by impact is 
reported in references 12 and 13. The SL'idy indicated that the influence of impact damage 
on ultimate strength may be negligible for lightly- lc .oderately-loaded panels but can be 
pronounced for heavily-loaded panels designed to carry loads at high strains. The present 
paper reports additional damage tolerance studies on compiession panels designed to carry 
compression loads without buckling and to match the cxtenfional and shear stiffnesses of 
typical commercial aircraft wing panels. The stiffness requirement results in designs 
which are heavier than those required to meet buckling requirements alone. The present 
experimental '^'•udies included lightly- and heavily-loaded blade-stiffened panels as well as 
heavily-loaded hat-stiffened panels. 

A photograph of the experimental apparatus is shown in fig^ire 19. Specimens were proof 
tested to the design load to ensure panel integrity. Some specimens were impacted with no 
applied axial load and then loaded in compression to obtain their residual strength, while 
other specimens were impacted while loaded. Some loaded panels that did not fail catas- 
trophically upon impact were subsequently loaded to failure to obtain their residual 
strength. 

Test parameters and results of this study are presented in Table V and a summary of 
the test data on the effect of low-velocity impact damage on the ultimate strain of stiff- 
ened panels loaded in compression is presented in thj bar grap'' of figure 20. Selected 
data taken from reference 13 are also shown in figure 20. The designs range from lightly- 
loaded (.53 MN/m) panels to heavily-loaded (5.8 MN/m) panels. The horizontal dashed lines 
in figure 20 represent the theoretical buckling strain at the indicated design load and « 
hat or blade symbol is above multiple test data of a corresponding design. The solid 
circles indicate the applied axial strain at the time of impact. If a panel failed on 
impact, the circle is at the top of the bar. If a panel did not fail on impact, the 
bar extends above the solid circle to the failure strain measured during the residual 
strength test. 

Blade-stiffened panels .- Strength tests on impact-damaged panels were conducted on 
blade-stiffened specimens constructed of designs D, E, and G (Table II). Lightly-loaded 
specimens Dl, El, and E2 exhibited postbuckling capability without failure during proof 
tests. Specimen Dl was impacted while under load at twj locations between stiffeners at 
velocities of 101.0 and 99.6 m/s without failure. The axial strain at impact for these 
two tests was 0.0015 and 0.0020, respectively. Extensive back surface damage was inflicted 
as can be seen in the photograph of panel Dl shown in figure 21(a). The damaged panel was 
capable of carrying a load in the postbuckled state that was approximately 28 percent 
greater than the undamaged panel buckling load. A photograph of the moire-fringe pattern 
illustrating the skin buckle pattern is shown in figure 21(b). Test results for specimens 
El and E2 indicate these lightly-loaded panels were also unaffected by these impact damage 
test conditions. Specimen E2 sustained complete penetration by the projectile and was 
capable of carrying a load 24 percent greater than the panel buckling load. These test 
results and results reported in references 12 and 13 indicate the ultimate strength of 
lightly-loaded panels designed to carry load at relatively low strains (less than 0.003) 
are essentially unaffected by low-velocity impart damage of the type studied. 

Results for heavily-loaded blade-stiffened specimens (design G) indicate the ultimate 
strength of these panels was degraded by impact at velocities around 100 m/s. During 
proof tests, design G panels were subjected to loads as high as 4.91 MN/m 
(strain 0.0068). Following impact at 61.0 m/s at zero load, panel G2 carried a load 
of 5.59 MN/m (strain ■ 0.0075) without failure. Further impact in the adjacent bay at 
92.4 m/s at zero load, however, caused the panel to fail at 3.94 MN/n axial load 
(i train ■ 0.0052) . 

Panel G4, impacted at 100.0 m/s while loaded to 2.62 MN/m (strain ■ 0.0035), had a 
residual strength of 2.79 MN/m which corresponds to an axial strain of 0.0037. The panel 
failed by delamination of the skin with stiffeners remaining undamaged. The series of 
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TABLE V.- THE EFFECT OF IMPACT-DAMAGE ON COMPkESSION PANEL STRENGTH 
(a) Blade-Stiffened Panels 
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Projectile used was 1 .27-cm-diametcr aluminum sphere. ‘^BuckleJ without failure. 

'Lette- indicates design described in Table II. Sa.or post buckling deformations at failure. 








(a) BACK SIDE DAMAGE 


(b) SKIN BUCKLE PATTERN 


1010 m/s I 99.6 n«A> 
STRA«N»O.OOI5 I STRAIKNa0020 


IMPACT SITES 


FIRST IMPACT SECOND IMPACT 


Figure 21.- Post buckling of a lightly-loaded 
blade-stiffened panel with impact damage. 
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Figure 22.- Propagation of impact damage 
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(o) BLADE- STIFFENED PANEL. (b) HAT-STIFFENED PANEL. 

Figure 23.- Ultimate failure of heavily-loaded panels damaged by impact. 

three moire-fringe photographs shown in figure 22 illustrates the progressive propagation 
of the failure boundary which occurred as the load approached Che ultimate value. The 
inset on the first photograph shows the C-scan of the damaged panel taken prior to the 
residual strength test. The C-scan indicated a damage region approximately 6.3 cm by 3.4 
cm. The failure of panels C2, G5 and G6, involved delamination of the stiffeners as well 
as Che skin. A photograph of panel G2 following failure is shown in figure 23(a). Based 
on these test results, design G impact-damaged specimens have a failure threshold 
strength and strain of about 2.79 MN/m and 0.0037, respectively. 

Hat-stiffened panels .- Impact-damage strength tests we-.e conducted on heavily-loaded 
specimens of designs H and I (Table II). No significant differences were observed in 
the failure load of impact-damaged panels of design H (all tape) and design I (mixed 
tape and fabric). Prior to impact testing, panels of both designs were proof loaded to an 
axial load of 4.05 MN/m and axial strains of 0.0055. For these hat-stiffened panels 
damaged by impact at approximately 97 m/s, the failure threshold strength and strain 
respectively, were greater chan 2.99 MN/m and 0.0042 (see Table V(b)). A photograph of 
panel H2 following failure is presented in figure 23(b). The failure initiated in Che 
impact region and involved delamination and shear failure of Che skin as well as destruc- 
tion of the stiffeners. 

Effect on Structural Efficiency 
of Limiting Ultimate Design Strain 

Pending improvements in the damage tolerance of composite structures and the develop- 
ment of better analytical failure prediction techniques, limiting the ultimate design 
strain appears to be Che most effective means for insuring the integrity of composite 
structures. Although test data are insufficient to establish an accurate limiting strain, 
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the value for the 4S0K cure graphite-epoxy material studied herein is likely to be less 
than 0.004. 

To establish the effect that limiting the design strain would have on the structural 
efficiency of graphite-epoxy composite wing panels, a study was conductec using the PASCO 
code. Results are presented in figure 24. The cross-hatched region represents the mass 
requirements for current aluminum commercial aircraft compression wing panels. The lov/er 
curve represents the mass requirements for a graphite-epoxy, blade configuration 
designed to meet the representative aluminum wing stiffnesses listed in Table 111. The 
graphite-epoxy designs also carry the axial compression loads without buckling assuming an 
initial bow-type imperfection (f- « +.003). Although the lower curve had no upper limit on 
strain, a value of 0.0093 was Che highest calculated strain required for the load index 
range studied. Curves are also presented for graphite-epoxy panels with these same con- 
straints except the maximum axial strain has been limited to either 0.003 or 0.004. 

For a value of of 0.7 MPa, a mass savings of nearly 50 percent is predicted for 

the graphite-epoxy design cotapared with Che aluminum design. For N^/L equal to 5.0 MPa, 
mass savings of approximately 45, 30, and 15 percent are predicted for graphite-epoxy 
designs with no limiting strain and maximum strain values of 0.004 and 0.003, respectively. 
Considering the complete wing compression cover, the total mass effect of limiting the 
design strain is an accumulative effect found by considering Che mars requirements for all 
Che different loads and regions of Che wing. Results for one such study are reported in 
reference 14 in which it was found chat limiting the graphite-epoxy design strain to values 
of 0.004 and 0.003 resulted in mass savings of 26 and 19 percent, respectively, compared 
to the aluminum designs. 


Damage Tolerance Improvement by Alternate Matrix 

The mechanisms of damage initiation and failure propagation usually involve matrix 
failure either through delamination or intraply matrix fracture. The implied matrix dependency 
suggests increased damage tolerance may be achieved by use of ai alternate matr: t material. A 
preliminary survey of commercially available epoxy preimpregnated materials identified 
one possible epoxy candidate. The effect of impact damage on the strength of this material 
was studied using plate specimens constructed of a 48-ply {±45/02/±45/02/±45/0/90)2g 
laminate. The laminate and test conditions permit direct comparison of this graphite-epoxy 
(material B) with Che graphite-epoxy (material A) data presented earlier, since the same 
graphite fiber was used in both materials. Lamina material properties for material B 
are listed in Table 1. Laminates of material B were 18 percent thicker than laminates 
of materin’ A, primarily due to the higher resin content (38 percent versus 28 percent). 

Results of this investigation are presented in the strain versus impact energy graph 
of figure 25. Data for material A is reproduced from figure 14. As before, solid 
symbols represent specimens loaded to tne indicated strain that failed oh impact, and open 
symbols represent the strain applied to specimens that did not fail on impact. The 
failure threshold strain is substantially higher for material B. For example, Che failure 
threshold strain for a 100 m/s impact speed is 0.0028 for material A compared with 0.0062 
for material B. 

It is important to establish the material properties 'hat affect laminate damage 
tolerance. Some insight may be obtained by studying the neat-resin tensile properties. 
Neat-resin tensile properties for materials A and B are shown in the inset of figure 25. 
Material \ has a slightly higher tensile modulus (4.0 GPa versus 3.1 GPa for 
material B) . The ultimate strength of material B is approximately twice Chat of 
material A and the strain at failure is approximately four times as great. Studies 
conducted on matrix materials with tensile stress-strain characteristics similar to those 
of material B indicate additional properties such as shear modulus and strength, resin 
content, and strain rate sensitivity may also be important factors affecting the damage 
tolerance of composite structures. 
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Figure 24. ~ Effect of limicing strain on structural efficiency 
of graphite-epoxy wing panels and comparison with aluminum wing panels. 
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Figure 25.- Damage tolerance improvement using alternate matrix. 
Projectile used was 1 . 27- cm-diameter aluminum sphere. 
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CONCLUDING REMARKS 


To achieve the maximum potential structural efficiency for composite panels, a thorough 
understanding of the factors that affect structural performance including buckling and 
material strength is required. Recent advances in analysis and sizing procedures reported 
herein permit composite panel designs to satisfy buckline requirements with greater 
accuracy and therefore improved structural efficiency. The effect of transverse shear 
deformations has been included to improve the accuracy of the buckling analysis for 
sandwich- blade-stiffened panels. A correction factor (based on infinitely wide orthotropic 
plate solutions) has been applied to improve the solution accuracy for low wave-number 
inplane shear buckling. Furthermore, a stiffened panel sizing code called PASCO has been 
developed that uses accurate analyses for buckling and an efficient optimization algorithm 
for rapid sizing. Additional analysis refinements allow PASCO to account for bow-type 
initial imperfections, lateral pressure, and thermal effects. The code is used to illus- 
trate the effects that selected configurations, loadings, and special stiffness constraints 
have on composite panel structural efficiency. 

Test results for stiffened panels designed by these improved procedures correlate well 
with theory when factors that affect buckling such as thickness and material property 
variations, initial imperfections, and transverse shear effects are adequately considered. 
Test results, however, also show that both low-velocity impact damage and circular holes 
can severely degrade the compressive strength of heavily-loaded composite panels designed 
to carry load at high strains. The problem is not critical for lightly- loaded designs 
since these panels usually carry the required load at relatively low strains. Low-velocity 
impact can cause local delamination and intraply cracking of the laminate. The local impact 
damage propagates in compression-loaded composite panels either by delamination, by local 
shear failure, or by a coupling between the applied axial load and the local lateral 
deformation due to impact. Currently, limiting the ultimate design strain based on 
experimental results appears to be one approach for addressing the damage tolerance 
problem. Although significant mass savings can still be achieved by imposing strains 
lower than are traditionally applied to metals, the need exists to improve the damage 
tolerance of composite structures. Test data indicate the impact-damage tolerance of 
graphite-epoxy laminates may be enhanced through improvements in the matrix material. 
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